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We propose the physical implementation of a Grover-like search problem by means of Frenkel 
exciton trapping at a shallow isotopic impurity against a background of competing mechanisms. 
The search culminating at the impurity molecule, designated the "winner" site, is marked by its 
enhanced interaction with acoustic phonons at low temperatures. The quantum search proceeds 
with the assistance of an Oracle-like exciton-phonon interaction which addresses only the impurity 
site, via the Dyson propagator within the Green's function formalism. The optimum parameters of 
a graph lattice with long-range intersite interactions required to trap the exciton in the fastest time 
are determined, and estimates of error rates for the naphthalene doped organic system are evaluated. 
We extend analysis of quantum search to a fluctuating long-range interacting cycle (LRIC) graph 
lattice system. 

PACS numbers: 03.65.Yz, 03.65.Ud, 03.67.Mii 

INTRODUCTION 

Quantum search via Grover's algorithm [l[ provides one possible framevifork by which combinatorial search can 
be speeded up considerably over those employing classical rules. For a function f{x) : {1,...,A^} — ?> {0,1} vifhere 
f{x) = for X ^ w and f{x) = 1 for x = w, the algorithm locates the "winner" w using of order ^/N queries for spatial 
dimensions exceeding two [2|. The successful bid is generally represented by an oracle with Uw\x) = (— 1)-''^'^^ |a;), 
and is sometimes treated as a black box which returns the value 1 for a particular one or a small set of inputs. The 
possibility of widening the choices for quantum searches via use of quantum random walks was first proposed by Farhi 
and Gutmann Q who showed that the time taken to move from one point of a graph vertex to another is exponentially 
faster for a walker following quantum rules over one who obeys classical physics. Unlike their classical counterparts, 
paths in quantum walks can interfere leading to nonexistence of a quantum walker at a vertex. 

In the continuous-time quantum walk formulation a Hamiltonian with oracle-like properties, = 1 — |ti')(i/7| 
is employed with jw) as the unknown state. Under this formalism, the ground state of the Hamiltonian is taken as 
\w) while all other equivalent but orthogonal states are considered as undesired outcomes of the search. Instead of 
the counting procedure, the efficiency of the search is quantified by the time taken to locate ground state 

In this work we attempt implementation of Grover's search via the Frenkel exciton, a bound system of electron 
and hole which propagates from one site to another in a quasi one-dimensional system of N molecular sites. Frenkel 
excitons are also modeled as delocalized excitation over the real crystal space [H| , and as noninteracting fermions , 
the exciton is thus an ideal example of an extended entangled system. We consider the exciton transfer mechanism 
as analogous to a quantum walker on a graph G where each vertex point represents a lattice site. The route of 
propagation of the exciton from one site to another is given by the edge joining one vertex to another with the 
number of edges incident on a specific vertex considerably being simplified for one dimensional systems. During 
a quantum search on graphs, a particular vertex is distinguished from all others and a successful outcome occurs 
when this vertex is reached. We examine the physical implementation of Grover's search by designating an impurity 
molecule as the "winner" located at the special vertex. The impurity site thus mimics the role of the desired outcome 
\w). The impurity molecule is distinguished from other molecules via the assistance of lattice vibrations, so that the 
difference between the energy band with and without the impurity present is accounted for by phonons. We note the 
fine balance involved in locating the impurity molecule, while lattice vibrations are needed to define vital properties 
of the impurity molecule, they should not be strong enough to hamper the search. The analysis can be extended to 
the existence of more than one impurity atoms, however we restrict our case to a single "winner" . 



EXCITON TRAPPING MECHANISM 

The trapping of excitons in molecular crystals has been well documented in the literature The overlap of 
intermolecular wavefunctions is weak in organic molecular crystal systems such as anthracene and naphthalene [sj, 
electronic wavefunctions are thus localized at individual molecular sites. The exact form of these functions are not 
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needed to derive salient properties of Frenkel excitons, however every molecular site has an equal probability of 
being excited by an incident photon due translation symmetry, a salient feature which has important implications for 
quantum parallelism. While a Frenkel exciton can be visualized as being smeared out in space, it can also be modeled 
as a bound electron-hole quasiparticle hopping from one site to another. 

Trapping sites are categorized as either impurity molecules (chemical traps) or dislocated molecules (physical traps), 
with the trap state assumed to be lower than the initial exciton state by an energy known as the trap depth. The 
excess energy is absorbed by phonons associated with lattice vibrations, the trapping mechanism thus occurring as 
a result of coupling between the exciton and phonons. A transition from a delocalized exciton state to a localized 
exciton state occurs at specific molecular sites. The extend by which neighboring lattice structure is altered depends 
on the magnitude of trap depth. A shallow isotopic impurity participates in the dynamics of the delocalized states 
of the host lattice system, unlike a deep impurity state. The energy state of the deep impurity falls far beyond the 
energy band states of the host molecular system. In this work, we consider introduction of a shallow isotopic impurity 
which alters the crystal symmetry only marginally, but one that eventually traps an exciton in the final state. We 
thus assume that the trap depth is present in the range of phonon bandwidths, so that a single phonon separates 
the host band and impurity trap. The trap depth may be viewed as a point of maximal lattice relaxation, with the 
excitation moving from an unrelaxed lattice structure to the relaxed lattice state at the impurity site, the surplus 
energy being converted into lattice energy. The discrete exciton spectrum in a finite one dimensional crystal system 
differs from the continuous spectrum of the infinite ID system. Hence a single discrete level corresponding to the 
self-trapped state remains distinct in the infinite chain unlike the discrete energy levels of the finite chain. 

We note that an initially delocalized exciton examines all of the sites simultaneously in a superposed form before 
reaching its target. The exciton lifetime therefore yields a measure of the efficiency of search from a quantum- 
information perspective. The absorption of a photon prior to search triggers the initialization process whereby an 
equal weighted superposition of all states is created. The exciton eigenvector appear as 

|K) =iV-i/2^e'K-'s/|0) (1) 
1 

where K is reciprocal lattice vector, |0) is the vacuum state with all molecules in ground state and is the creation 
operator of exciton at position coordinate 1. N is the number of unit cells in the crystal which increases with the 
number of molecules for each unit cell. We note that Eq. ([T]) is the result of a single step of a naturally occurring 
process involving photon absorption that require the equivalent of log2(iV) -I- 1 qubits to store the superposed state. 
In Eq. ([T]), we have dropped notations associated with the exciton spin, however we will consider spin indices in 
mechanisms involving exciton formation and annihilation in forthcoming sections of this work. 

The Frenkel exciton creation operator with wavevector K can be easily obtained using the Fourier transform of Eq. 
© 

Bl = 7V-i/2^e'K.i^t (2) 
1 

The exciton creation operator localized in K-space is delocalized in real space . The motion of the exciton in 
molecular systems is governed by an Hamiltonian derived using a tight-binding model 0, 0, II] 



E 



+ V A,m bIBi + V Mi,mBlB^ (3) 



where AE, the on-site (intra-site) excitation energy at equilibrium is the same at all sites due to translational 
symmetry. Di m is the dispersive interaction matrix element which determines the energy difference between a pair of 
excited electron and hole at a molecular site and ground state electrons at neighboring sites 0, [3|- Afi,m the electron 
transfer matrix element between molecular sites at 1 and m. Using Eq. ([2]), we obtain a Hamiltonian diagonal in K 
space 

He. = Y.Eo{k)BlBi, (4) 

k 

£;o(k) = AE+"^ Do^m + Mo,mexp{ik.m) (5) 
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where i?o(k) is energy of the exciton of wavevector k and subscript denotes the absence of lattice site fluctuations. 
For simple d— dimensional lattice with nearest neighbor transfer energy Mi „i ~ J, we obtain the well known dispersion 
relation Eo(K.) = AE + 2J^'^^^ cos{Ki) with the exciton band extending from AE — B to AE + B. B = 2Jd denotes 
the half band width and we emphasize its critical dependence on efficiency of exciton transport on designated graph 
structure of molecular sites. 



Excitonic bandwidth in graph lattice with long-range intersite interactions 

The bandwidth B which is determined by the energy span of the host lattice over a spectrum of wavevectors, 
is a vital property of the exciton. It influences excitonic coherence quantities such as its delocalization length and 
relaxation mechanism, and provides a measure of time needed for the exciton to move away from its initial site of 
excitation. It is intimately linked with arrangement of superpositions of single-site excitations of the initial state. In 
a recent work on photosynthetic organic systems Q , it was shown that the efficiency of excitation transfer depends 
both on the symmetry properties of such superposition states and on the number of sites among which the excitation 
is initially delocalized, with results showing an optimal delocalization length for which transfer time is a minimum. 
To illustrate the dependence of the trapping time on size-effects associated with long-range excitonic coherences as 
well as underlying structure of the lattice points, we consider transfer matrix element Mi^rn of the form 

Ml,™ = (6) 

\m — 

where Mm,m = and intersite energy term J > 0. This form of the matrix element has recently been shown to reveal 
various rich and interesting properties for a select range In the absence of lattice vibrations, 

the eigenstates of Eq. ^ form waves with momenta K = within the Brillouin zone fc € [— tI '^i^h energies 

ri=l ' ' 

where Li^(a) is the Polylogarithm function and ^'(z,/x,a) is the Lerch transcendent function [lo| . In the limit 
where K which occurs at the upper band edge, we use a simplified dispersion relation obtained as Ek ~ 
Eo{/J.) — JA{fi)\K\'^~^ where Eo{/j) is the upper band-edge energy and A{^) is a dimensionless constant. £'o(m) 
and A{^) can be obtained using series expansion of Eq. ([7|). An approximate but simple relation between the 
bandwidth and system size N follows at the upper band edge (1 < ^ < |) 

B{K -^0,N^oo)^ Eo{pi){l - AifiYN^'") (8) 

where ^(/i)' < 1 is a small dimensionless quantity. While the states are extended over the whole crystal structure at 
/Lt < I, it has been shown that localization effects appear for higher values of fi for systems with uncorrelated diagonal 



disorder and long-range coupling of the form in Eq. (|6]) [11[. We note that lattice vibrations, as considered in our 



work, provide a fabric of correlated diagonal disorder and therefore localization effects may be present with properties 
different from those obtained by Rodriguez et al [ll|. Nonetheless, we consider that the system obeying long-range 
coupling of the form in Eq. ([B]) still undergoes vital changes at the critical point when /i = | 

ESTIMATING SEARCH TIME USING GREEN'S FUNCTIONS. 

Implementation of the quantum search process involves the following crystal Hamiltonian with the impurity at site 



w 



Ht = Ho + H^ (9) 

Ho = Hex +^fiuj{q)blbq 
q 

H^ = ApBlB^ (10) 
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where Hex is given in Eq. Q and the second term in Hq denotes the phonon energies. Hq correspond to a impurity- 
free lattice structure. denotes a trapped exciton at site w, with trap depth Ap defined as the difference between 
a state of zero phonon and one where interactions between exciton and phonons leading to absorption of a phonon 
is optimum. We consider the presence of just one impurity site and assume that the translational symmetry of the 
molecular system is left intact. This is justified for isotopic impurities which mimic the molecular structure of the 
host molecules quite well. We note that the lowered energy Ap associated with the Hamiltonian H^, effectively marks 



the searched state, As discussed in Ref. [ij], solving the Grover's problem involves the determination of the 

lowered energy, and hence energy of the Hamiltonian in Eq. ^ . This translates to evaluation of the trapping time 
of an exciton initially located at a site other than at \w). 

We emphasize that the driving term fia;(q)6jj6q associated with a background of boson bath assist with the 
Oracle search operation. A resonant coupling term associated with exciton-phonon interaction at all host sites has 
been intentionally left out in Eq. (j9l). This is based on the assumed presence of an Oracle-like mechanism which 
ensures that all phonon-assisted relaxation (one-phonon emission and Raman processes) cumulates at the impurity 
site, bypassing other lower energy sites, during the search process. An Oracle is a component of an algorithm which 
can be regarded as a " black box" that recognizes solutions to a given problem [l^ . The inner dynamics of the Oracle 
is generally unknown, however at very low temperatures its inner workings appears almost apparent. Raising the 
temperature of the crystal medium therefore involves the risk of the Oracle malfunctioning, as will be shown in detailed 
evaluations in the next Section. We therefore address the impurity trapping problem at very low temperatures (close 
to OK), a limit in which the mean- free path associated with phonon scattering is large due to a linear scaling of 
exciton scattering rate with temperature. The exciton is more likely to be trapped at the impurity site than to bypass 
it. Various other competing mechanisms which act to disrupt the search and introduce errors will also be discussed 
subsequent to the current Section. 

The exciton-phonon coupling Hamiltonian associated with the interaction between the exciton at site w and phonons 
can be written as ^ 

i?e.p=i?iB.„5]x('Z)(&^ + foq) (11) 

q 

where coupling function x{<i) arises due to acoustic phonons dominant at low temperatures. x(g) is dependent only 
on the phonon wavevector q 

X{qf = ^hLo(ci)ELR (12) 



Elr — 
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where phonon frequency w(q) = ''^I'zIj ^ being the velocity of sound in the crystal medium and E^n is the lattice 
relaxation energy. / is the mass coefficient of host molecules, taken as almost the same as that of the impurity and 
deformation potential Eu remains an invariant throughout all lattice sites. A coupling coefficient 7 = Elu/B yields 
a measure of the scattering regime of the crystal system. For strong lattice vibrations at high 7, the exciton is less 
likely to hop and tends to remain localized at its initial site. At small 7, the exciton is in a delocalized state, until 
it reaches the site at \w). Thus the fine interplay between E^n and B plays a crucial role in the efficiency of the 
quantum search. While lattice vibrations are essential in locating the impurity site, these should not be too strong 
to disrupt the intermolecular cooperativity which determines the search process. 

We emphasize two critical operations of the Grover's search which is satisfied by Eqs. (fTO)) and ([TT|). While Eq. 
([T0| provides a local operation which distinguishes the impurity site from other sites, Eq. ([Tl]) provides an oracle- 
like operation which addresses only the impurity site via the evolution mechanism present in the Green's function 
formalism. The Green's function for an exciton at time t is given by 0, [l6| 

^(k,^) = -i(0:nq|T{Sk(t)[/(t,0)Bj,(0)}|0:nq) (13) 
Bk(0 = Sk(0)exp(i£;(k) t) (14) 

where |0 : riq) denotes a state with riq phonons of wavevector q and zero exciton population. T is the time ordering 
operator and [/(<, 0) is the Dyson function 

[/(t,0) = Texp [^i exp{iHTt')Hexpe^p{-iHTt')dt'^ (15) 
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where Ht and H^xp are given in Eqs. ([9]) and dTTI) . A direct relation between the Green's function G'(fc, t) in Eq. 
and exciton energy -E(k) has been derived by Craig et al flTj at t > 0, 

G(k,t) = -iexp[-i £;(k) (i^)] (16) 

At temperature close to T = OK, {F) can be obtained in an explicit form [l3, [3] 



q 



1 , , , 
h ^—(1 - exp(iryk,qt) 



(17) 



?7k,q = fta;(q) + (k)-Ap (18) 

is the mean energy of the exciton band, N is the number of sites and the pre-factor is due to the search for 
just one winner. 

The imaginary component of the Green's function in Eq. (IT6l) yields damping rate of the exciton transfer process 
among lattice sites. Here we use the damping rate to estimate the time needed to search for the impurity molecule, 
Ts 

^ = ^EW'?)I''^('?m) (19) 
* q 

In order to obtain an explicit result for T,, we assume that the phonon modes obey the Debye distribution with 
density of states given by p{huj) = ?>N {huj)'^ / {fujj dY' where lod is the Debye cutoff frequency. Using Eqs. ([8]), (fT2)) 
and (|19p. we obtain an approximate expression for search time at 1 < < | 

11 1 



(20) 



To 



AW 



N^'-' (22) 



where Eq{^) is the upper band-edge energy in Eq. (jS]). In order to highlight the explicit dependence of Ts on N , we 
have approximated energy of a single phonon, fvjj{c\) « Ap instead of using the energy conservation relation in Eq. 
(1181) . The search time constitutes a A'^-independent term associated with the maximum bandwidth and dependent on 
the topological parameter jj, (see Eq. ([5])), and a second Grover-like N dependent term which scales as TV"^. Thus two 
concurrent mode of search mechanisms appears to be involved in Eq. (j20l) , one associated with a quantum coherence 
effect that is independent of A^ and another term which has the well-known Grover-like search features [l[. At A^ — ?> oo, 
we obtain a constant search time Ts = Tq for the system. We note that the maximum bandwidth associated with 
Eo{ij) in Eq. ([5]) corresponds to a delocalized state of optimum coherence when the search is conducted from all 
lattice sites. At the peak of the energy band, the exciton can be considered to exist at all sites at the same time, in 
a maximally entangled state of existence. The search time Tq associated with the maximum bandwidth yields the 
fastest time and consequently the most efficient route to search the impurity site. For very large Tn > To, Ts ^ Tq, 
the search route is dominated by the first coherent term in Eq. ([20l) . 

At the marginal = |, we obtain a term that equivalent to the well-known Grover's search term, cx N2 [ij. The 
Grover's algorithm result gives the maximal possible probability of finding the desired element [l3|- However it only 
describes a specific mode to search a list of A^ sites and in this regard, the speedup features in the excess search term 
T/v of Eq. pn)l at /i < I is associated with the peculiarities of the level spacing (of order Ai? ^ A^^+^) at the top 



band as discussed in Ref. [ll|. This implies that the search proceeds with a decreasing number of target sites as 
the energy of the quasi-particle decreases. The explicit relation in Eq. ()20|) was obtained by utilizing a background 
boson field coupled with the evolution mechanism of the Dyson propagator (see Eq. (fT5|)) within the Green's function 
formalism. Such an approach may facilitate search procedures with a number of channels for the delocalized excitation 
to pursue the preferred impurity site efficiently. The range of ji considered here provide ideal conditions required to 
trap the exciton in the fastest time, when the quantum search is considered accomplished. In contrast to the dynamics 
at the top of the exciton band, it was shown that the energy level spacing scales as /S.E ~ A^^ at the bottom band 
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This means that the search for a target site located at a higher level would experience a slowdown effect as is 
obvious when substituting higher ^ values in Eq. The parallels between the quantum search mechanism and 

physical processes considered so far is shown in Table 1. 
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Table 1 Parallels between exciton trapping mechanism and the proposed Grover-like search problem 
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Quantum Zeno Effect 


Error Minimization 



Organic molecular systems with /i = 3. 

It is important to note that the marginal exponent /i = | is not expected in any known molecular system, nev- 
ertheless in organic systems such as Frenkel excitons in Langmuir-Blodgett films, molecules experience dipole-type 
exchange mechanisms with /i = 3 [^, 0|- For such systems, the exciton bandwidth expression in Eq. ([5]) is modified 
as fl2| 

B{K)^Eoi^l){l~A{^lyN-^) (23) 
for > 3). Using Eqs. ([23]), ^ and (dU, we obtain 

. (25) 

where Tq is given in Eq. (|20p . We use crystal parameters for the naphthalene host crystals (CiqHs) doped with 
isotopic impurities used in an earlier work [20|, with Ap « 50cm~^, huo ^ ^o(m) ~ 90cm~"'^ and Elr — 0.004 eV. 
For N = 10^ and /i = 3, we obtain search times Tq ^ O.Olps, ~ 10 — lOOps and overall Tg « O.Olps. The vast 
difference in search times, Tq and T/v suggests that the search route in organic systems may be dominated by the 



first coherent term in Eq. ([Ml) . Delayed fluorescence techniques [2l| may be used to detect the impurity site and 



confirm the results obtained here. The emission of fluorescence during experimental measurement can be considered 



to signal the end of the quantum search. Recent advances in spectroscopic techniques [2^ may be used to verify the 
femtosecond time scale estimates predicted for Tq in naphthalene doped organic systems. Studies of exceptionally 
efficient excitation transfer in photosynthetic systems has also been explored in recent works 0. In this regard, 
there is possible application of the coherent size-independent, connectivity dependent time component Tq derived in 
Eq. ([20]) to studies related to excitation transfer in photosynthetic systems and associated phenomenon of quantum 
photosynthesis. 



PROCESSES WHICH COMPETE WITH GROVER-LIKE SEARCH PROCESS 



In this section, we consider various sources of decoherence-type mechanisms that will result in improper search 
patterns of the impurity site. 
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Downward relaxation and decoherence of superposed exciton states 



Exciton-phonon scattering can interfere with the search process by causing downward relaxation of states which 
occupy the upper exciton band levels. The possibility of such a process can be minimized by ensuring that the exciton 
remains delocalized and hence entangled over a number of molecular sites with size equivalence to the overall number 
of sites to be searched. In the presence of decoherence effects due to lattice vibrations, it is thus imperative that the 
coherent length associated with the superposed exciton state in Eq. ^ should at least exceed the number of sites to 
be searched. This means that the exciton state should retain its coherent form at least for the time needed for the 
exciton to reach the target impurity site, a condition best described by the relation 

Tn < r,cat (26) 

where Tgcat is the scattering time. Using perturbation theories, the scattering rate can be estimated using Tgcat ^ 
Elr^bT p{E) where p{E) is the normalized density of states of the unperturbed lattice system[83. For s constant 
density within the band {IS.E — B, AE + B), p{E) = jj^2_p.2 ^ — i- The density of states, p{E) is highly 

sensitive to the parameter p which appears in Eq. (jG)) 12"]. For p = 3, we obtain 

hB 

where E^j^ the lattice relaxation energy is defined in Eq. and B denotes the half band width. Comparing Eq. 
(p7)) with Eqs. (|2T|) and ([22]) . and substituting the upper band-edge energy Eo{p) = 2B and A{py w 1, we obtain 
an explicit condition required for the Grover-like search to proceed successfully without being affected greatly by 
downward relaxation 

The number of sites that can be searched without disruption increases at low temperatures and at higher trap depth, 
Ap. We bear in mind, however the possibility that the impurity energy band may become disengaged from the exciton 
band with significant loss in translational symmetry due to deep traps associated with very large Ap. 

Another undesirable effect associated with exciton-phonon scattering is the introduction of composite exciton- 
phonon characteristics to the exciton bandstructure. It can be seen that a fixed sum of exciton and phonon wavevectors 
given by K = ^ (k -f- q), involves an infinite number of phonons with wavevector q. This leads to a weakening 
effect of an impurity trap due to interference. For instance, the presence of an isotopic impurity displaces the energy 
band by £' = B^ + Ap^ , so that at bandwidth B — 300cm~^ and Ap « lOOcm^^ we obtain an energy shift of about 
AE = 16cm~^ below the band-edge. At temperatures higher than 30K, it is not uncommon for appearance of band 
spectra features of the order ~ 10cm~^ Q, this implies that exciton-phonon coupling leads to comparable energy 
shifts as the impurity state. Therefore exciton-phonon relaxation may produce "defective" winners. Hence operating 
at the right temperatures holds the key to maintaining a fine balance between undesirable features of phonons and 
its facilitating role in energy transport. The calculations above that at temperatures close to 5K, we expect a clear 
distinction between "defective" and "genuine" winners when the exciton are trapped at the impurity states displaced 
from the exciton band. Another side effect of Tgcat is manifested in the the phase coherence in Eq. ^ which becomes 
increasingly damped at higher temperatures 

Bl = A^-i/2^e^(K+ik').i5t ^ ^-i/2^g^K.ig^kM^t (29) 
1 1 

It is well known that the imaginary component fc', which appears as a result of the scattering mechanism in Eq. (j27p . 
is a source of decoherence in the entangled exciton state. 

For the naphthalene host crystals (CioHg), typical value of i? ^ lOOcm^^, and using ksT = 3.7meV at T = 30K, 
we obtain Tgcat ~ 50ps which reduces to 500ps at a lower T — 5K. These times are comparable to the search time 
T/v ~ lOOps obtained earlier, suggesting that downward relaxation mechanism may pose as a serious drawback to the 
Grover's search mechanism at increasing temperatures (T > lOK). It should noted that Tq (see Eq. pT|) ) associated 
with the maximum bandwidth and with coherence properties still remains the dominant pathway to energy transfer 
due to the very fast times To ~ O.Olps computed earlier for naphthalene crystal system in comparison to the Grover's 
search term T/y (see Eq. ([121)) as well as Tgcat- 
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Exciton formation at impurity site 



It is possible that an exciton is created at an impurity site due to the binding of free charge carriers instead of 
trapping of an exciton participating in the Grover-hke search mechanism. In this section, we compute the rates 
associated with such a process which acts to disrupt the search process. In Frenkel excitons, both charge carriers are 
generally created via optical excitation on the same molecule. The transition from the initial state of charge carriers 
to a final state of trapped exciton via acoustic phonon occurs as 

e 

e + h — > I + huj, (30) 
h 

e 

where e, h and ^ denote a free electron, hole and exciton respectively, and haj is the frequency of an acoustic phonon. 
The initial and final states appear as 

N> = E4(^e)4-k(^/«)|0;^> (31) 

l/> = Bl_^{<js)bl\0;n) 

where ae{(Th) denotes the electron (hole) spin, while the lowest internal exciton state is assumed with spin as = 0(1) for 
the singlet (triplet) state [1^. ci\j^d\^ denote the electron (hole) creation operator of wavevector k in their respective 
energy bands. The exciton wavevector can be written as a linear combination of all eigenvectors of the free charge 



carriers (see Eq. 3 in Ref.[23|). The initial and final state eigenvalues are obtained using 



F -W , ^'k^ , ^!(K_lk)! , p .o2^ 

_ W,+E, + ^^+ 2rn* + ^^'^ ^^^^ 

where Tph represents the phonon energy that remains unchanged during the transition and Wa is the total electronic 
energy of crystal with all molecule in their ground states. m*,m'^ and M* are the respective effective masses of the 
electron, hole and exciton. The rate of formation of the exciton at the impurity site is evaluated using the following 
exciton-impurity-phonon interaction operator 

-,1/2 

(k - k') • e(q) exp[z(k' - k + q) • p] 4Sk(6^ + feq), (33) 



k,k',q 



2IpU}{q) 



where e is the unit polarization wavevector and Ip is the mass coefficient of the impurity molecule which we assume 
to the same as the mass coefficient of the host molecules. The rate of formation of an exciton at the impurity site is 
calculated using 

k,q 

where ni denotes the concentration of free charge carriers. For the singlet exciton {as — 0), we obtain using Eq. 
([3T|) ■ ([32]) and (|34|) an explicit expression for the exciton formation rate at low temperatures 

+ (35) 

where N* and do are the number of molecular sites and inter-molecular distance respectively at an edge of the crystal 
unit cell. Et is the transport gap, the minimum energy required for the formation of a free electron and hole pair. 
El, is the binding energy, ^ is mass density of the host molecules and v denotes the velocity of sound in the crystal. 

For naphthalene host crystals (CioHg) doped with isotopic impurities, we use Ap « 50cm^^, {Et ~ Eh) = 3 x 
10-*cm-\ n, ~ lOi^cm-3, ^ = 1283 kg m^^, y ^ IQ^cm s''^, N x da ^ 7.8 x 10^1° m, and obtain Tf « lO'^s. 
We note the crucial dependence of the formation rate on the carrier concentration, however even with a large carrier 
concentration n^, a comparatively slow error time is obtained. This can be attributed to the large bandgap in molecular 
systems, with the requirement that charge carriers have to be optically generated the probability of exciton formation 
at the impurity site is less likely at low excitation intensities. Moreover due to the large binding energies of the Frenkel 
exciton, participation of phonons of relatively higher energies is needed before excitons can be formed at the impurity 



site, dissimilar to the situation in semiconductors [24|. We can therefore conclude that exciton formation at impurity 



site is not disruptive to the proposed Grover's like search mechanism, at least in naphthalene -doped organic crystals. 
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Exciton-exciton annihilation. 



Excitons in close proximity can dissociate due to electrostatic interactions, creating free charge carriers. The energy 
released from recombination of an electron and a hole bounded together as an exciton at one molecular site can lead 
to dissociation of exciton at another site, with possible participation of phonons. Such a process can compete with the 
proposed Grover-like search process, hence we examine exciton-exciton annihilation between two free excitons prior to 
exciton trapping at the impurity site. For simplicity, we consider the two interacting excitons to be in the singlet spin 
state ((75 = 0), the generated hole carrier to be trapped at a site leaving only the election free after annihilation. The 
conservation of momentum rule ensures that the exciton momentum is transferred to electron after the annihilation 
process, in which participation of phonons are omitted. Due to the small phonon participation number, the order of 
magnitude of the annihilation rates are not greatly affected by omission of interaction terms associated with phonons. 
We write the eigenvectors of the initial and final states as 

N) = Bl^BljO:n) (36) 
l/> = E«k(^'=)4'K)|0;^) (37) 

where (^^(rfk) which denote the electron (hole) creation operators are associated with free charge carriers. The 
transition from \i) to |/) is enabled by electrostatic interactions Vmn between molecules m and n with a transition 
matrix of the form 

{f\V\i} = 2iV-i^exp[k.(n-l)] (38) 

l,n 

e e e e 

(0,m;e + h,l |Kri„| I , n; I m) - -(e + h,l;0,n |Kn„| I , n; I , m) 
h h 2 h h 

The first term within the square bracket describes recombination of a free exciton (denoted by ) at site m, with 
the released energy inducing the exciton at n into a pair of free carriers (denoted by e+h) at site 1. A similar 
interpretation holds for the second term, with recombination of the free exciton taking place at site n instead. The 
initial and final state energies are obtained as 

= Wo + 2EoiK) + ^^ (39) 
Ej = Wo + E,{K) + ^ (40) 

Wo is the total electronic energy of crystal with all molecule in their ground states and Eo{K) is the exciton energy 
given in Eq. ([S]). Ee{K) denotes the free electron energy. Using the Golden Rule formula of Eq. (IM)) as well as Eqs. 
and p9p . we obtain the exciton-exciton annihilation rate Ran as 



1 10 n. 



T hi 

an If' 



=^3/2 7-^5/2 j\r=^6 7 6 
TO ' EJ N do 



(41) 



where n^x is the concentration of free excitons and N* and do are defined below Eq. p5|) . Ex, the transport gap is 
used to estimate the square of the transition matrix in Eq. psp. This approximation is justified as Eq. psp describes 
transition associated with the energy required for the formation of a free electron and hole pair as quantified by the 
transport gap. 

We estimate the exciton-exciton annihilation rate Ran for the naphthalene host crystals (CioHg) using Et = 3 x 
lO^'cm-i, Hex ^ lO^^cm-s^ ^ ^ 1283 kg m'^, iV x do = 7.8 x 10'^° m, and obtain Ta„ = ~ IQ-^s. Thus exciton- 
exciton annihilation time Tan is comparable to the Grover-like search time, Tjv ~ lOOps at high concentration of free 
excitons. At even higher exciton concentration n^x > 10^'^cm^^, this process may dominate the exciton dynamics 
leading to a breakdown of the Grover-like search process. Hence it is imperative that the exciton concentration is set 
low so that exciton-exciton annihilation is not disruptive to the search process, in naphthalene-doped organic crystals. 
We note the crucial dependence of Ran on N* , which yields a measure of packing density of the searched sites. Hence 
the distribution of searched sites plays an important role in the influence of exciton-exciton annihilation, which is 
expected intuitively as well. 
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Exciton fusion and fission. 



During exciton fusion, two triplet excitons combine to form a singlet exciton as shown below 




The subscripts T and S correspond to the triplet and singlet exciton respectively, with excess energy generally 
released as molecular vibrations. In naphthalene organic system, the energy difference between two triplet excitons 
and a singlet exciton is given by 2E{as = 1) — E{{as = 0) = Ey— 1.3 eV, with the triplet exciton energy tending to 
decrease faster than the singlet exciton .25i |. A reverse process in which a singlet exciton undergoes conversion into 
a pair of triplet excitons is termed exciton fission. Hence in the fusion (fission) process, a singlet (triplet) exciton 
is formed with all excess energy is released (absorbed) via intra- or inter- molecular vibrations. The calculation of 
rates of exciton fission or fusion is similar to that for exciton-exciton annihilation for which free charge carriers were 
created in the final state (see Eq. (HH])). We therefore do not give quantitative details of the process, however we 
note the sensitive dependence of such processes on riga;, the concentration of free excitons as obtained in Eq. (j4ip . 
This was also noted in a recent work with rates estimated in the picosecond range ~ 1 — 10 ps for closely related 
naphthalene-related systems [2^. We expect that at low exciton concentration and packing density of the searched 
sites N* , exciton fusion and fission processes are less likely to disrupt the overall dynamics involved in the Grover's 
like search process. 



EXCITON DYNAMICS CONTROL VIA THE ZENO EFFECT. 



Quantitative estimates for the naphthalene doped organic systems show that several mechanisms, in particular 
downward relaxation due to exciton-phonon interaction, exciton-exciton annihilation, exciton fusion and fission process 
become dominant at raised temperatures and high exciton concentration which translates as packing density of the 
searched sites. Overall we note that the coherent term Tq in Eq. pT|) associated with the maximum bandwidth 
still remains the dominant pathway to energy transfer in comparison to the Grover's search term T/v (Eq- (|22p ). 
relaxation term Tgcat (Eq. p7p ). formation time Tf (Eq. (j35|) ) and exciton-exciton annihilation time. Tan (Eq. 
(j4T|) ). Nevertheless, the exciton dynamics associated with Grover's search term T/v can be carefully harnessed via 
the quantum Zeno effect (27j . This effect embodies the controlling effect of a measurement process whereby a system 
monitored to determine whether it remains in a particular state resists making transitions to alternate states. In recent 
years, this idea has been studied using an adiabatic theorem ^2^^ '29'] in which different outcomes are eliminated and 
the system evolves as a group of exclusive quantum Zeno subspaces within the total Hilbert space. The quantum Zeno 
effect can be utilized to exclusively direct the Grover-like search process to proceed uninterrupted while suppressing 
alternate mechanisms of exciton transport and dynamics. Such directed control of many-body dynamics via Zeno effect 
has been experimentally demonstrated in recent works (sol . Isij , with potential applications in quantum information 
systems. We are unaware of the use of Zeno effect to switch selected interactions "on and off' in molecular systems. 



SEARCH ON A FLUCTUATING LONG-RANGE INTERACTING CYCLE (LRIC) GRAPH LATTICE 

A long-range interacting cycle (LRIC) graph is an extension of the nearest-neighbor model and a simpler variant 
of the model considered earlier in Eq. ([B]) . In a LRIC graph lattice, all two nodes of distance m are connected [ll] , 
a key feature which ensures that all lattice sites have the same connectivity of 4. These newly added edges serve as 
shortcuts in the cycle graph, and serve to increase the cross paths and hence bandwidth B of the lattice system. The 
one-dimensional LRIC graph lattice is denoted by two parameters, N and m and the exciton Hamiltonian present in 
such a lattice appear as 

Hex = +'^Y. +'^Y. ^/^'-l + 5Z B\Bi+,n + JJ2 BjBi^rn (43) 

I I I I I 

where Ei denotes the onsite energy and J > 0, and the intersite energy is assumed to be independent of the distance 
between any two sites. Using the Fourier transform in Eq. ([2]), the eigenenergies are obtained as En^m = Ai? -|- 
2J{cos K + cos mK) with momenta K = within the Brillouin zone fee [1,^]- The LRIC Hamiltonian is diagonal 
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in K space 

i?e. = (44) 

K 

The bandwidth _B of a Frenkel exciton in a LRIC lattice therefore has a simple form i? = 2 J cos + 2 cos miiT. We now 
consider that the lattice points of a LRIC graph undergo small vibrations, and using the Green's function approach 
considered in the earlier Section, we can examine the effect of m on the search time. We consider a large N so that 
the impurity trap falls within the bandwidth B. The exciton bandwidth is also considered to include the acoustic 
phonon bandwidth. At small m, B{K ^ 0, iV -J> cxd) « 2J(1 - (^f) and using Ea.([T9|. we obtain a search time 
component Tjy that scales as N'^ . However for large m, the bandwidth oscillates strongly depending on m. Using 
iV = pm where integer p provides a measure of the lattice network extension, and considering small values of cosmic, 
we obtain 



To 



cos(f) 



where as in Eg. (1201) . Tq is an TV- independent term. The search time component Tp decreases gradually with p so 
that search time increases overall as is expected with an increasing lattice network extension. We note that 
Ts in Eq. (|45l) does not resemble the Grover-like form obtain earlier in Eq. (fT9| . due to the vast reduction in 
exciton delocalization property associated with the simpler exciton Hamiltonian structure of the LRIC graph lattice. 
Nevertheless the explicit dependence of search time Tg on p illustrates the importance of the efficiency of excitation 
transfer in quantum searches. The search time Ts is optimum for large m, which agrees well with the result that the 
exciton survival probability decays faster on networks of large m (32| . 

DISCUSSION AND CONCLUSION 

We have demonstrated the plausibility of a Grover-like search problem by means of Frenkel exciton trapping process 
at a shallow isotopic impurity. We have considered several mechanisms which may compete on equal terms with 
our proposed Grover's search mechanism. Quantitative estimates for the naphthalene doped organic systems show 
that Grover's search mechanism remains plausible under suitable conditions and may be physically implemented with 
assistance from Zeno effect. We emphasize the striking resemblance of a derived explicit term for exciton trapping time 
to the well-known Grover's search time, which includes notable speedup features in the characteristic range 1 < /i < | 
for a graph lattice with long-range intersite interactions. These parameters provide optimum conditions required to 
trap the exciton in the fastest time. We suggest use of Zeno effect associated with experimental measurements as 
a viable means of eliminating ubiquitous disruptive interactions and harnessing the search process, with emission of 
fluorescence signifying the end of a quantum search. 

A recent work on light-harvesting networks [s^ showed the role of network size and connectivity as critical factors 
which lead to optimum energy transfer in graph lattice systems, and our results agree well with this observation. We 
emphasize the significance of a derived connectivity dependent time component To, associated with intermolecular 
cooperativity and coherence properties of the entangled exciton as a dominant source of excitation transfer mechanism. 
The derived explicit term can be utilized to examine high efficiencies of energy transfer in photosynthetic systems 
. Estimates of various search times and excitonic mechanisms obtained in this work can be verified by experiments 



employing fast femtosecond time resolution spectroscopic techniques [22|. Finally, our results show that the natural 
evolution mechanism of the Dyson propagator within the Green's function formalism may play a critical role in the 
construction of efficient quantum search devices, and in bridging the gap between excitonic theory and its use in 
quantum information processing. 
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